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Abstract 

We investigate the influence of higher twist corrections to deep inelastic structure functions in 
the low-Q^ and small-x HERA region. We review the general features of the lowest-order QCD 
diagrams which contribute to twist-4 at small-x, in particular the sign structure of longitudinal 
and transverse structure functions which offers the possibility of strong cancellations in F2. 
For a numerical analysis we perform a twist analysis of the saturation model which has been 
very successful both in describing the structure function and the DIS diffractive cross section at 
HERA. As the main conclusion, twist 4 corrections are not small in Fl or Ft but in F2 = Fl+Ft 
they almost cancel. This indicates the limitation on the use of the DGLAP formalism at small 
X and Q^. We point out that F^ analysis needs a large twist-4 correction. We also indicate the 
region of validity of the twist expansion. 
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1 Introduction 



A deeper understanding of the transition from perturbative QCD to nonperturbative Pomeron 
physics in deep inelastic scattering at low and small x remains one of the central tasks in HERA 
physics. Approaching the transition region from the perturbative side, one expects to see the on- 
set of large perturbative corrections - in particular those which belong to higher twist operators 
in QCD. The twist expansion defines a systematic approach and, therefore, provides an attrac- 
tive framework of investigating the region of validity of the leading-twist NLO DGLAP evolution 
equations. The essentials of the theory of higher twist operators and their Q^-evolution have been 
laid down twenty years ago: a choice of a complete operator basis has to be made [Q], and for the 
evolution [|| one needs to compute evolution kernels which, for partonic operators in leading order, 
reduce to 2 — > 2 kernels. The problems of mixing between different operators of a given twist has 
also been addressed in |^. Explicit calculations have been done mainly for fermionic operators. 
In the small-x region at HERA, however, we expect gluonic operators to be the most important 
ones. Recently, a first attempt has been carried out to analyze the twist-4 gluonic operators in the 
double- logarithmic approximation (DLA) H],^. In addition to analytic calculations also a first 
numerical analysis has been presented. As one of the main results, it has been pointed out that, 
due to a complicated sign structure, subtle cancellations among different twist-4 corrections are 
possible. As to the numerical results, the freedom in choosing initial conditions for twist-4 gluonic 
operators, in combination with our presently very limited knowledge of the twist-4 evolution equa- 
tions, make a systematic QCD study of higher twist corrections in the low-Q^ and small- a; region 
at HERA a rather difficult but challenging task. 

In order to gain a detailed insight into the role of higher twist it may be helpful to discuss in some 
detail the simplest low-order QCD-diagrams (rather than using the whole Q^-evolution machinery 
collected in [^). Particular attention has to be given to the question of possible cancellations 
between different contributions. As to the choice of input distributions, it seems advisable to make 
use of more specific, model dependent assumptions on the input distribution. The most reasonable 
starting point, in our opinion, is that model which has been most successful in describing the low-Q^ 
data of HERA: the saturation model of |^ which contains only four free parameters. This model 
not only describes very well the 7*p-cross section in the low transition region where the role of 
higher twist is of particular importance, but also allows to connect, in a quantitative way, the total 
cross section data with the DIS inclusive diffractive process A particular benefit of using this 
model is the interpretation in terms of QCD diagrams: comparing with the analysis of the QCD 
diagrams it is possible to read off a choice of twist-four initial conditions. Since the model (before 
doing any twist expansion) describes the HERA data, it is also likely to provide a realistic estimate 
of twist-4 contribution in the low-Q^ and small-x region. 

We begin with reviewing the expressions for the simplest QCD diagrams and discussing patterns 
of possible cancellations. In the following part we review the saturation model and define the twist 
expansion. In the third part we perform a numerical analysis and draw our conclusions on the the 
magnitude of gluonic twist-4 corrections. The results of our analysis are in qualitative agreement 
with the estimates presented in , suggesting that twist-4 corrections to F2 = Ft + are small 
down to Q2 ^ 1 GeV^, X 10 ^. However, we also find that this smallness is due to an almost 
complete cancellation of the twist-4 corrections to Ft and Ft: both of them, individually, are large, 
but have opposite signs and nearly the same magnitude. We interpret this as a warning against 
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(a) 



(b) 



(c) 



FIG. 1: The simplest QCD diagrams with 2,3 or 4 t-channel gluons. 



using the twist-2 formahsm at too low and small x. 



2 QCD Diagrams 

We are interested in twist-4 corrections AFj to the transverse (T) and longitudinal (L) structure 
functions 



(2.1) 



A QCD analysis of twist-4 corrections at small-x starts from the lowest order diagrams shown in 
Fig. 1. The photon can have transverse or longitudinal polarization. Approximate expressions for 



the fermion loop, -D, 



{T,L)(abcd) 
4:0 



can be found in they are valid for small x, and corrections are 



of the order 0{x). The diagrams with two gluons start with leading twist, but they also contain 
higher twist - similar to the BFKL approximation which also contains corrections of all orders in 
Unfortunately, we have no way to obtain information on the higher- twist couplings to the 
proton. The diagrams with three gluons, through the reggeization of the gluon, are higher order 
corrections to the diagram with two gluons; in the analysis of [Q] they are needed to complete the 
covariant derivative of the two gluon diagram. In a complete twist-4 analysis these diagrams with 
two and three t-channel gluons have to be included, but presently we do not know how to estimate 
their magnitude. The most interesting twist four diagrams, presumably, are the ones with four 
gluons. They belong to the four gluon operator which is expected to play the most crucial role at 
low and small x. Its coupling to the proton has been discussed in |^, where arguments have 
been given that the simplest model which respects the AGK cutting rules consists of (at least) two 
pieces: 



abed 
iP4 



abed I abed 
¥'43 + f4A 



(2.2) 



where 



, abed 



3-8 



r\j-^ A/^ 



d'^'d'^fsil, 2; 3, 4; u;) + <5'^'=<5'"*/s(l, 3; 2, 4; u) + 5'^'*r/s(l, 4; 2, 3; co) ) (2.3) 



jac xbd 



ad xbe , 
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and 



^fK = - — pppri (r''"rtA(l,2;3,4;.;) + 

+ r^"^r''/A(l,3;2,4;..) + r'^-/-''^/A(l,4;2,3;a;)) . (2.4) 

Here /s and /a are w-dependent, positive- valued functions which play the role of unintegrated 
gluon densities. Combining them with the quark loop expressions, 

AF^, = -TT^ Df^"^' (^)-=^ , (2.5) 



and retaining only those terms which give rise to Q^-logarithms, we arrive at the following twist-4 
corrections to the transverse and longitudinal structure functions: 

and| 

(I)) - ^-*<-)i ■ <-) 



The w-dependence of initial conditions (p4s{^) and ipiA{^) will be assumed to lead to a power-like 
behavior of the form (1/x)^^ where the exponent A is unknown. Thus, together with the l/Q^ 

We changed the definition of Fl in comparison to Q . Now, Fl is twice the previous one to have F2 = Ft + Fl 
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suppression, these twist-4 corrections are of the form 

Q 



From this general observation one immediately sees that the value of where twist four becomes 
important is j;-dependent. One of the most striking features is the sign structure: the transverse 
and longitudinal cross sections, (2.6) and (|2.7]), have opposite signs, and in AF2 = AFt + AFl one 



faces a strong cancellation. If 9943 and 994A are of the same order of magnitude (such that the square 
bracket expression is positive), we expect the twist-4 corrections to F2 being slightly dominated by 
the negative corrections to Fl, i.e. the higher twist corrections to F2 could be small and negative. 

It is important to note that these corrections to the deep inelastic structure functions are closely 
related to the twist-4 corrections to the cross section of diffractive qq production (the s-discontinuity 
line between gluon 2 and 3). Twist-4 in diffractive qq production has been observed experimentally, 
and it has to be a part of the twist-4 corrections to Ft or Fl. However, contrary to the most naive 
expectation, the two-gluon systems on both sides of this cutting line cannot be restricted to be in 
color singlet states: the AGK rules in perturbative QCD require a structure of the form (2.2) 
and (2.3), and from this one easily sees that the system of gluons 1 and 2 is not restricted to color 
singlet. 

Turning to corrections of the form a^lnQ^/Qo ^^S- ^^^^ * mixing problem. In |^] it 

has been argued, on the basis of an all-order analysis in the double-logarithmic approximation, 
that there are several gluonic twist-4 operators, four-gluon operators and one two gluon operators. 
They have different Q^-evolution equations, and they couple to the proton with different couplings. 
Furthermore, they mix, i.e. there are transitions from two-gluon states in the t-channel to four- 
gluon states. When computing aglnQ"^ /Qq corrections to the diagrams in Fig. 1(c), we group the 
contributions in the form shown in Fig. 2: (a) illustrates the transition of the two-gluon state to 
the four-gluon state and leads to the correction AF^, and (b) illustrates the first evolution step 
of the two-gluon operator and gives the correction to Ai^^. The calculation of these diagrams is 
described in and leads to the following twist-4 contributions: 
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Note by the comparison of (|2lll), (pl^) with (^), (13), respectively, that AF^ gets additional 



QglnQ^/Qo corrections without changing the structure of initial conditions. The corrections AF^ 
are of the order a| (05 In Q^/Qq)) thus they are not present in the lowest order result which is 
proportional to a|. 

For low Q^-values it is not a priori clear whether these corrections to the twist-4 contributions 
are important or not: there is an additional suppression factor and for low Q^-values the 

logarithm log Q'^/Qq does not provides much enhancement. To get a first idea, it may, again, be 
useful to draw a connection with diffractive dissociation. As illustrated in Fig. 2(a), these diagrams 
describe diffractive production of qqg systems. There is no doubt that these diffractive states have 
been observed at HERA: a direct analysis of their twist-4 component (e.g. the observation of 
diffractive final states with only hard jets) would provide a direct evidence for the presence of these 
higher twist corrections in the deep inelastic structure function. 

A simple analysis of twist-4 corrections could be based upon the presented low-order expressions. 
However, even within this framework we need two initial conditions, cps and (/)a- Relating them to 
the twist-4 diffractive qq cross section (as described in Q) gives only one condition, and, hence, is 
not enough. We are therefore lead to build a model for the initial conditions. The most successful 
description of the low-Q^ transition region at HERA has been provided by the saturation model 
of B, and we will use this model to determine the initial conditions. 



3 Twist Four in the Saturation Model 

Let us first briefly review the model of and its decomposition into twist components. It is well 
known that the 7*p-cross sections, 

^t,l(x,Q2) = l^FT,L(x,g2), (3.1) 
can be written at small X ctS loL IcJ '. 

aT,L{x,Q'^) = j d^r dz\^T,L{^,r)f a{x,r'^) (3.2) 

where ^T,Liz,'r) denotes the transverse and longitudinally polarized photon wave functions, and 
a{x,r'^) is the dipole cross section which describes the interaction of the qq pair with the proton. 
In addition, z is the momentum fraction of the photon carried by the quark, and r is the relative 
transverse separation between the quarks. The wave functions are solely determined by the coupling 
of the photon to the qq pair, see e.g. In the dipole cross section is assumed to depend on x 
through the ratio of the transverse separation r and the saturation radius Ro{x), and the following 
form is proposed: 

a{xy) = -o<7(^) - -o{l--p(-^)} • (3.3) 

At small r (r <^ 2i?o), the dipole cross section grows quadratically with r, cj ~ aor'^/'iRQ, while for 
large r (r ^ 2i?o), it saturates, a = (Tq. In order to describe the energy dependence both of the 
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total DIS cross sections and of the low-Q^ cross section measured at HERA, the saturation radius 
has the following x-dependent form: 



with Qq = 1 GeV. Eqs. (|3.2| )- (|3.4| ) define the saturation model. The physical motivation for such a 
parameterization and its significance for diffractive processes in DIS is discussed at length in [^,0. 
The three parameters in the model are determined from a fit to the total DIS cross section data 
at X < 0.01 and look as follows: do = 23 mb, xq = 3 • 10"'* and A = 0.29. In this way a very good 
description of data in a broad range of and x is obtained. In fact there is a fourth parameter in 
the model, an effective quark mass ruf = 140 MeV in the photon wave function, chosen such that 
the results of the model, extended down to photoproduction region, are in a good agreement with 
photoproduction data measured at HERA. 

In order to evaluate the cross section (|3.2| ) it is convenient to employ the Mellin transform to 
factorize the wave function from the dipole cross section: 



°° dv f 1 \ V2+i!^ / rnl 



In the case of zero quark mass we obtain 

u( n^ _ 6«em ^ 2 ^ 9/4 + ^ ^ ^ sinh(^z.) r(3/2 + iv) 

HT{i^^^)-^L.^fTQ^T^\T^M^)) vr^ _r(-i/2-.^.) ^^-^^ 



and 



rr I _ 6aem 2 ^ 1/4 + TT ^ sinh(7r;^) r(3/2 + it/) 



/ 

In addition, Giu) in ( |3.5| ) equals: 

~2 /'"2N-3/2-ii^ ,^2 



G{y) = / df'^ (f2) g{f^) = -r(-l/2 - ii/), (3.8) 

Jo 

for the saturation form of the dipole cross section g{r^) = 1 — e^^'^. From relation ( |3.5| ) we see that 
the cross sections ax l depend on x in the saturation model only through the combination 



{3S 



^ Q^Rlix) Q^\xJ • 
Notice that a similar combination, eq. (|2.8D, occurs in the twist-4 analysis in the previous section. 
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FIG. 3: The solid line ~ 1) in the (x, Q^) plane indicates a critical line of the model The area between 
the two dotted lines corresponds to the acceptance region of HERA. 



The detailed discussion in |Q (which will not be repeated here) shows that there is an essential 
change in the energy (and also Q^) dependence of the cross section if we move from the "perturba- 
tive" region, ^ < 1 (large and not too small x), to the nonperturbative "Pomeron" region, ^ > 1 
(small and very small x). In the former region we have a power-like rise in 1/x (in agreement 
with the observed rise of the structure function at small x), whereas in the latter one the cross 
section stays constant (apart from a logarithmic enhancement factor). Hence the region ~ 1 (see 
Fig. 3) marks the transition from one region to the other. In the region ^ < 1 it is natural to 
expand the inclusive cross sections in powers of ^ (which means powers of while in the region 

^ > 1 in powers of or Q^. It is clear, however, that the expansion in powers of 1/Q^ cannot be 
valid down to = 0. Therefore, approaching the transition region ~ 1 from the perturbative 
side one expects the leading-^ approximation to fail somewhat to the right of the line ^ = 1 (called 
a critical line in [^). 

With expressions (|3.5| )- (|3.9| ) the above discussion can be made more precise. The expansion in 
powers of ^ will be identified as "twist expansion" (below we will explain why). The power series 
in ^ (or 1/^) is determined by the singularities of the integrand in eq. (|3.5| ) in the complex z/-plane. 
These are single or multiple poles, located ai v = zbz(2n + l)/2 with n = 0, 1, 2, .... 

The z^-integration runs along the real axis, and for < 1 it is tempting to close the contour in the 
lower half-plane. A closer look reveals that in such a case an asymptotic expansion for our cross 
sections is constructed. This is done by computing residues of the poles in the lower half-plane 
which leads to an expansion in powers of ^ or 1/Q^. Now, the critical line (or better, a strip) 
indicates a limit on validity of the asymptotic twist expansion. It is also possible to close the 
contour in the upper half-plane in which case an convergent expansion in positive powers of 1/^ or 
is obtained for any value of ^ 7^ 0, 00^. However, this expansion is not practical in the large 

^We thank Lech Mankiewicz for a discussion on this point. 
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or small ^ analysis. 



The first singularity encountered in the lower plane is a pole at u = —i/2. In the saturation 
model the transverse cross section has a double pole which generates a logarithmic behavior for the 
leading-twist contribution: 

CTT = ao ^ H-C - mi + e ln(l/0) , (3-10) 

The longitudinal leading-twist contribution has only a single pole and therefore does not produce 
a logarithm: 

a^ = ao^e^^e (3.11) 



Higher-twist contributions can be obtained by evaluating the residues at the lower lying poles , the 
twist-4 contribution at the pole u = — 3i/2, etc. The results for the transverse contributions are: 

Twist-4: 

/ 

Twist-6: 

- = -1:4^ (ilje' - ime + lenvo) . (3.13) 

Twist-8: 

— E 4^ + |^(^)«' - . (3.14) 

For the longitudinal contributions we find: 
Twist-4: 

a, = a^Y^ef-^ (-|^^ + ^^(3^^ - ^^^ln(l/o) , (3.15) 

Twist-6: 

— i:e}^(S«'->^'-|«'>H' 

Twist-8: 
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It is not surprising to see a strong similarity between the twist-4 contributions ( p. 61) and 
computed to lowest order in QCD, and those found in the saturation model ( ^.121) and ( |3.15|) . In 
particular, the sign structure is the same: cjt is positive while 0"^ is negative. Also, contains a 
logarithm while ax does not. Finally, the ratio of the twist-4 contributions, 

94/75 + 4/5 (ln(l/0 - V;(3)) 
^ 375 ' ^^-^^^ 

is similar to the ratio obtained from the lowest order QCD calculation 

AFf _ 94/75 + 4/5 ln(QVQ2) 



AF^ 3/5 



(3.19) 



Thus, the saturation model can be interpreted as a result of summing the diagrams shown in 
Fig. 1(a) and 1(c) (and more "iterations of gluon ladders"): the leading twist comes entirely from 
Fig. 1(a), twist four from Fig. 1(c) etc. Viewed in this way, the model can be used to define values 
of the initial conditions for ( |2.6| ) and (2.7), but only in the combination 14(/94g — 9(/94a. The higher 



order corrections illustrated in Fig. 2(a) are not included in the model. The success of the model 
might indicate that this type of corrections is not very important near the transition (critical) line; 
but it may also be that a suitable modification of the model might be necessary. 



4 Numerical Analysis 

Further insight will be gained through a numerical study. We choose to use the expressions of section 



3 and compare the full expressions (3.5) (with m/ = 0), which can be integrated numerically, with 



the twist-expansion ( 3.10D -( ^.17| ). As far as the twist-4 terms alone are concerned, we could have 



started from the saturation model, extract the initial conditions and than turned to the QCD 
approximations listed in section 2. For our present discussion, however, we find it more instructive 
to study the role of higher twist (in particular, twist-4) in a more general context. The saturation 
model, which describes the data, provides an analytical formula for the cross sections and allows 
to investigate the twist expansion and its breakdown near ^ = 1. Since the dependence upon 
and X is through the variable ^, we present the numerical results as a function of ^. 

In Fig. ^ we show lines of constant ^(^ = 2, 1, 0.1): variation of ^ means moving from one line to 
another. With the help of Fig. 3 ^ can be translated into the x and variables. We begin with 
presenting the ratios of the leading-twist approximations to the full cross sections as a function 
of ^. In Fig. 4 we show the ratios for the transverse a'Jf^'^/ax, longitudinal a'^j^'^/uL, and the 
total (o"£^^ + aTf^'^)l{aT + cl), cross sections. The striking result is that for both the transverse 
and longitudinal cross sections separately the higher twist corrections become large when we are 
approaching the transition line ^ = 1. Moreover, the transverse higher twist correction is positive 
while the longitudinal one is negative. In the sum, however, those corrections almost cancel each 
other and the overall correction is small. 

This is a possible explanation of the smallness of the higher twist correction to F2 found in the 
analysis of MRST |l^. The authors of this analysis use the following simple parameterization: 

■?HT I ^ r<l\ T?LTr r^2\ ( ^ , ^2{x) 



{x, Q') = (x, Q^) ( 1 + ) , (4.1) 
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X 


D2{x){GeV'') 


0-0.0005 


0.0147 


0.0005-0.005 


0.0217 


0.005-0.01 


-0.0299 


0.01-0.06 


-0.0382 


0.06-0.1 


-0.0335 


0.1-0.2 


-0.121 


0.2-0.3 


-0.190 


0.3-0.4 


-0.242 


0.4-0.5 


-0.141 



TABLE 1: The MRST values of the higher-twist coefficient D2{x) in Eq. (4.1) 



and determined the function D2{x) from a fit to DIS data. The result is given in Table 1 which we 



reproduce from |1C]. In general, for small x, D2{x) is small and negative but it becomes positive 
for the smallest x. We found a similar result, the leading twist approximation deviates from the 
exact formula by less than 10%. The sign structure of this deviation also agrees with the MRST 
analysis; it is negative but to the left of the transition line the deviation becomes positive. The 
last result should be taken with some care since, strictly speaking, the twist expansion for ^ > 1 
in the saturation model makes no longer sense, and a new expansion in powers is appropriate. 
We extrapolated, however, the leading twist formula to that region and found an agreement with 
the exact result up to ^ = 2. This indicates that the phenomenological success of a leading- twist 
analysis might be deceptive: the leading-twist approximation to F2 remains a good approximation 
also in the region in which the whole twist expansion has already collapsed. 

The last point is also illustrated in Fig. 5 where the individual higher twist components are shown. 
We plot (T^^ -|-(t£ for r = 2, 4, 6, 8, together with the exact result (solid line). The overall impression 
is that near ^ = 1 all higher twist corrections are getting large, leading to the conclusion that the 
concept of higher twist becomes meaningless. If we naively extrapolate the higher twist formulas to 
the region ^ > 1 they diverge. Nevertheless, to the right of = 1 there is a region where twist 6 and 
8 are small and can be neglected, whereas twist-4 accounts for the deviation between the twist-2 
approximation and the exact result. Fig. 6 illustrates this in another way: the sum of twist-2 and 
twist-4 provides a rather accurate description of the exact formula. In this region (in Fig. 3 between 
the two lines ^ = 0.2 and ^ = 0.9) twist-4 corrections should improve the QCD description of deep 
inelastic scattering. 

In Fig. 7 we present the results of the analysis performed only for the longitudinal twist contributions 
o"^ and r = 2, 4, 6, 8. The striking result in comparison to the total cross section analysis is a large 
and negative twist-4 contribution which accounts for a large difference between the exact and 
leading twist result. This is shown in Fig. 8: again the sum of twist-2 and twist-4 provides an 
accurate description of the exact formula. Notice that the difference — cr|^'*^* is large already 
at ^ = 0.1. Taking this result, we conclude that an analysis of the longitudinal structure functions 
Fl based entirely on the leading twist result is unreliable already for quite high values of Q'^ and 
not to small x. A similar effect, although not so pronounced, occurs for the transverse cross section 



(large twist-4 corrections for the transverse case were also discussed in [11]). The role of higher 



twist in the longitudinal structure function was broadly studied in |14] and a qualitatively similar 
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results to ours were obtained there for small x. 

In summary, the main lesson to be learned from this study is the cancellation of higher twist in F2: 
twist-4 is not small, neither in the transverse nor in the longitudinal part, but it is hardly visible 
in the sum of both because of the mutual cancellation. Moreover, the twist-2 approximation to F2 
works even beyond ^ = 1 where the whole concept of a twist expansion should make no sense. As 
far as is concerned, the leading twist significantly exceeds the exact result for ^ > 0.1, and a 
large and negative twist-4 correction is necessary to obtain an agreement with the exact result. 

5 Discussion and Conclusions 

In this article we have carried out a simple numerical analysis of gluonic twist-4 corrections in the 
low-Q^, small-x HERA. We have reviewed what an analysis of lowest-order QCD diagrams suggests. 
One of the most striking features are differences in sign between transverse and longitudinal twist-4 
corrections which may lead to a small twist-4 correction, even if the corrections to Ft or Fl are not 
small at all. Because of the unknown initial conditions more input is needed. We then have used 
the saturation model which provides an excellent description of both F2 and the DIS diffractive 
cross section at HERA, and we have analyzed its higher twist content. We found a one-to-one 
correspondence between the twist expansion of this model and the QCD diagrams discussed before. 
The model can therefore be used to define initial conditions of a QCD higher twist analysis. 

In our numerical analysis we have restricted ourselves to a careful study of the saturation model. 
We found that, indeed, twist-4 corrections to F2 remain small, and this smallness is due to an 
almost complete cancellation between large corrections to Ft and F^. This implies that although 
twist four corrections are small in F2, the use of the leading- twist DGLAP formalism for extracting 
structure functions becomes doubtful in the low-Q^, small-x HERA region. Within the saturation 
model we have quantified the limit of applicability. 

This problem is even more acute for Ft alone, where the large twist-4 correction is directly ex- 
posed, providing a crucial contribution which brings the leading twist result close to the exact one. 
The DGLAP formalism might not be reliable in this case even for higher values x and Q^, see Fig. 8. 

Clearly, our numerical conclusions are based upon a specific model. The phenomenological success 
of this model provides some reasons to believe that the conclusions are realistic. Moreover, on a 
qualitative level our conclusions are in agreement with the independent fit to the HERA data of 
| |10| . Nevertheless, some uncertainty remains. We have outlined that in the saturation model some 
contributions are not present which one would expect to see when starting from QCD diagrams. If 
included they may modify the subtle balance between transverse and longitudinal structure func- 
tion. They may also shift the transition region (the "transition line" of the saturation model, in 
reality, may turn out to be a rather narrow "transition strip"). However, the coincidence with the 
MRST fit makes us feel that the conclusions of our analysis are "not far from reality" . 
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transverse 




FIG. 4: The ratio of the leading-twist and the exact cross section in the saturation model : for the transverse, 
a'^'^/cjT, longitudinal, a}^'^ /ul, and the total, a^^j^/ cjt+l, cross sections. 
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total 




FIG. 5: The exact total cross section (solid line) and the higher twist contributions in the saturation model 
as a function of the parameter ^ defined in the text. The cross sections are rescaled by a common 
factor. 
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FIG. 7: The exact longitudinal cross section (solid line) and the higher twist longitudinal contributions in 
the saturation model as a function of the parameter ^. The cross sections are rescaled by a common 
factor. 
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FIG. 8: The same as in Fig. 7 but with the leading twist and twist-4 added (dotted line). 
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